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G ■ Abstract 

It is shown how the bihncar differential equations satisfied by Fredholm determinants 
r~^ I of integral operators appearing as spectral distribution functions for random matrices 

may be deduced from the associated systems of nonautonomous Hamiltonian equations 
satisfied by auxiliary canonical phase space variables introduced by Tracy and Widom. 
The essential step is to recast the latter as isomonodromic deformation equations for 
families of rational covariant derivative operators on the Riemann sphere and interpret 
(^ • the Fredholm determinants as isomonodromic r- functions. 

O" 

^ . . . 

O ■ 1 Differential equations for Fredholm determinants 

0\ ; in random matrices 

In the theory of random matrices, it is known that in suitably defined double scaling limits 
the generating functions for spectral distributions are given by Fredholm determinants 
of certain integral operators ^J El El El- For example, in the universality class of 
^ [ the Gaussian Unitary Ensemble (GUE), in the bulk of the spectrum, the probability 

of having exactly {mi, . . . ,mn} scaled eigenvalues in the sequence of disjoint intervals 
^. {([01,02],... , [a2„_i,a2n]} is 

5— J ' ( — 1)"* Qrh^S ^^ 

^^"^^' • • • '"^") = ^,!...rnJ9Ar---W" A,=...=A.=i' ^ = 2:-^' (1-1) 

3 

where r is the Fredholm determinant 

T^ := det (l - k^\ (1.2) 

of the integral operator Kg : L^(]R, C) -^ L^(M, C) with the sine kernel 
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Rescaling at the (soft) edge of the spectrum, the corresponding quantity is given by the 
Fredholm determinant 

r^ := det (l - K^^ (1.4) 

of the operator with the Airy kernel J18j 

Kv\{x) = 2^Xj -— v{y)dy, (1.5) 

j=i •^'»2j-i -^ y 

where Ai{x) is the Airy function. If the measure is taken to be the one associated with 
either the Laguerre or Jacobi orthogonal polynomials, rescaling at the (hard) edge leads 
to the Fredholm determinant 

r„^ := det (l - K^) (1.6) 

of the operator with Bessel kernel [HI El 

where Ja{x) is the Bessel function with index a. 

It was shown by Tracy and Widom |171 1181 [TU] , extending earlier results of the Kyoto 
school |llj . that all these Fredholm determinants can be computed by quadratures in terms 
of solutions of certain associated nonautonomous Hamiltonian systems in which the end 
points {flj} play the role of multi-time deformation variables. Moreover, these Fredholm 
determinants may be interpreted as isomonodromic r-functions 1161 IIUI [3] in the sense 

of [nmni. 

More recently, Adler, Shiota and van Moerbeke [21 [31 have shown that the Fredholm 
determinants r , r^ satisfy hierarchies of bilinear differential equations with respect to 
the endpoint parameters. These follow from combining Virasoro constraints satisfied by 
certain associated KP r-functions with the bilinear equations they also satisfy with repect 
to the KP flow parameters {ti,t2) • • •}) evaluated at the zero values of these parameters. 
The approach of [21 [HI was based on the application of vertex operators, integrated over 
the intervals {[a2j-i,a2j\}^ to suitable "vacuum" KP r-functions, effecting thereby a con- 
tinuous version of Darboux transformations, yielding new KP r-functions, such that the 
Fredholm determinant equals the ratio of the two. 

For the Airy kernel, the first equation in this hierarchy may be expressed as 

-D^pA _ /^j)j)^pA ^ 2Po-F^ + 6 [VlF^f = 0, (1.8) 

where 

F^:=lnr"^ (1.9) 



and 

2n 

da 



^--E^rTT"' "'^^' (^-^^^ 



j=i 



532 J Harnad 

while for the Bessel kernel, it is 

VfF^ - 2V\F^ + (1 - a^) vIfJ^ + Pa^iFf 

- \v2F!^ - 4 (PiFf ) {vIf!^) + 6 [VlFJ^f = 0, (1.11) 



with 



F^:=lnT^. (1.12) 



No analogous equations were derived for the sine kernel, although in the special case where 
the intervals [a2j-i, cl2j] are chosen symmetrically about the origin, the Fredholm determi- 
nant T'^ may be expressed ^JE] as a product rfr^i of two Bessel kernel determinants. 

2 2 

In the case of a single interval, it is easy to see that equations H1.8|) and ()1.1H) just give 
the T- function form of the Painleve equations Pji and Py, respectively, to which the Tracy- 
Widom systems reduce in the case of the Airy and Bessel kernels. It seems reasonable to 
expect that analogous results hold for the general case, involving an arbitrary number of 
intervals. The purpose of this work is to show how the hierarchies of equations derived 
in [21 El can in fact be deduced directly from the Tracy-Widom Hamiltonian systems 
for both the Airy and Bessel cases, and to also apply this approach to the sine kernel 
case. The main step is to recognize that the Hamiltonian systems imply isomonodromic 
deformation equations for associated families of rational covariant derivative operators on 
the Riemann sphere. It is known jl21 I13j that such isomonodromic deformations give rise 
to bilinear equations for indexed sets of isomonodromic r-functions related by Schlesinger 
transformations. The fact that for the systems associated with the Airy and Bessel kernels 
such equations may be written in terms of a single scalar r-function is due to the presence 
of a pair of conserved quantities, allowing the elimination of the additional variables by 
fixing the level sets of these invariants. In the sine kernel case this is not possible, and 
the associated bilinear equations therefore involve coupled systems for r together with 
a pair of additional variables {rf, t^). 

In Section 2, equations (|1.8|) and (|1.1H) are first derived directly from the Hamiltonian 
systems of |18U19j . In Section 3, it is shown how the isomonodromic deformation equations 
following from the associated Hamiltonian systems may be used to derive the full hierarchy 
of T-function equations for all these cases. In section 4, these results are related to the 
rational classical i?-matrix approach to isomonodromic and isospectral systems developed 
in PIHl. 

2 Deduction of T-function equations 
from the Hamiltonian systems 

To establish notation, following ^3^1^], we define the quantities: 

X2j := 2i^{I - i^)-V(«2j), 2;2j+i := 2yA;(I - i^)-V(«2j+i), (2.1a) 

2/2i := i\/A;(I - Ky'ijia2j), y2j+i := ^/\j{^ - K)" V(a2i+i), (2.1b) 



ra2, 



xo:=2^Aj/ 0(x)(I-K)^V(a;)dx, (2.1c) 
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2/0 := Yl ^i / '^(^^^^ - Kr'cP{x)dx, (2.1d) 

where, for the case of the sine kernel K = K , 

, , , sinfvrx) , . x / x , ^ 

(/,(x) := ^ — ^, iIj{x) := cos(7rx). (2.2) 

vr 

while for the Airy kernel K = K , 

cP{x):=Aiix), ,l;(x):=^^, (2.3) 

and for the Bessel kernel K = K^ , 

^{x):=M^), ^{x):=x^^^. (2.4) 

(An odd number of variables may also occur if we set one of the Oj's equal to some 
fixed constant, say or oo, and eliminate the corresponding pair {qj,Pj)-) As shown in 
J17l 1181 I19j , the logarithmic derivatives of the associated Fredholm determinants are given 
by: 

^s.^^0F^^^^2 ^lyi^m^yj^ (2.5) 

A J •'A ^—^ n- — n. 



duj 4 -^ -^ 4 f-^ Qj — au 



for the sine kernel, 

J fe=i J " 

for the Airy kernel, and 

ajG^j := aj-^ = ^i ~ 16 ^" ~ "^' ^ ^^-^ ^^' 

for the Bessel kernel. 

For use in what follows, we also define the quantities 

2n 

R^ := VmF^ = Y^ afG^, m e N (2.8) 

i=i 

for the sine kernel case, 

2n 

Ri:=V^F'' = Y,ci^Gf, men (2.9) 
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for the Airy case and 

2n 

<^ := VraFl^ = J2 afGlv ^^^ (2-10) 

i=i 

for the Bessel case. For aU three cases, we define the following sequence of bilinear forms 

2n 2n 2n 

P^-=Y.<yl Qm:=J2<^l S^:=Y,aJxjyj, m € N. (2.11) 

j=i j=i j=i 

As explained below, the {G^}'s and {G^ j's may be viewed as sets of Poisson commu- 
ting, nonautonomous Hamiltonians on an auxiliary phase space with canonical coordinates 
{xoiVoiXj^yj}, such that the quantities defined in (|2.1|) satisfy the corresponding systems 
of Hamiltonian equations. These equations will then be shown to imply equations H1.8|) 
and (P~TT|) . 

2.1 The Airy kernel system 

The system of dynamical equations for this case is given _18^ by 

3 + K (2.12a) 

J + k, (2.12b) 

+ 2yj-yoXj, (2.12c) 



+ 2^^j - ^o)xj + yoXjyj, (2.12d) 



dxj 


1 [xjyk - 


VjXk)xk 


dak 


2 ttj - 


-ak 


9yj 


1 {xjyk - 


VjXk)Vk 


dak 


2 a,- 


- ak 


dxj 


_ 1 v- (xjVk 


- yjXk)xk 


daj 


2^ a 

fc=i 


■j - ak 


dyj 


1 ^ {xjyk 


- yjXk)yk 


daj 


2^ a 


j-ak 


dxQ 
daj 


= -XjVj, 


dyo _ 1 

daj 4 



(2.12e) 



Viewing the a/s as multi-time parameters, and the quantities {xo,yQ,Xj,yj} as canoni- 
cal coordinates, this is a compatible system of nonautonomous Hamiltonian equations 
generated by the Poisson commuting Hamiltonians {Gf} defined in ()2.6() . There is an 
additional functionally independent Hamiltonian, defined by 

G^ := yl - xo - ^Qo, (2.13) 

which also Poisson commutes with all the G^s. Since Gq is not explicitly dependent on 
the parameters {oj}, it follows that it is a conserved quantity. Since all the quantities 
{xo,yo,Xj,yj} defined in ()2.1|) vanish in the limit {oj — > oo, V j}, the invariant Gq must 
vanish on this particular solution. Therefore we may express xq in terms of the other 
variables as 

xo = yl- \Qo- (2.14) 
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The quantity Rq defined in (|2.9j) will just be denoted 

R:=R^ = Y, Gf = Po - -Qi + -ylQo - yoSo - -QI (2.f 5) 

i=i 

where ()2.14() has been used. In terms of R, equation H1.8|l becomes 

V^R - iViR + 2R + 6{VoRf = 0. (2.16) 

It follows from the Poisson commutativity of the Hamiltonians {Gj}j=i^,,,^2n that their 
Hamiltonian vector fields applied as derivations to R give zero, and hence along any 
integral surface of eqs. 1)2. 12() . the derivatives of R with respect to the Oj's are just given 
by its explicit dependence on these parameters. This just comes from the Qi term in 
expression (|2.15() . and therefore we have 



aR __i 2 



-x2 (2.17) 



Comparing with (|2.12e|) . this implies that 

Gi:=yo-R (2.18) 

is a second conserved quantity. Since in the limit {uj -^ oo, V j}, both yo ^-nd R vanish, 
G^ must vanish for all values of the parameters, and therefore the invariant relation 

yo = R (2.19) 



is satisfied by this solution. Applying 
that 


the operators 


Vo, Vi to R, 


it follows from (|2.1*7() 


VoR = 


-\Qo, 












(2.20a) 


ViR = 


-\^^- 












(2.20b) 


Eqs. (|2.12I) also imolv that aoDlication 


ofPo 


to {Qo, 


So,xo,yo 


,Qi} 


gives 




^oQo = 


-- 45o - 2yoQo, ^o^o = 


2^1 


--X0Q0 + 2P0, 






(2.21a) 


VqXo = 


-So, Voyo = -^Qoi 












(2.21b) 


VoQi = 


= Qo + 4Si-2yoQi. 












(2.21c) 



Further application of Vo and Vi , using H2.20a() , 1)2. 21|) and (|2.14|1 , therefore gives 

V^R = ^yoQo - So, (2.22a) 

ViR = -igi + 2yoSo - ^ylQo - \qI - 2Po. (2.22b) 

Substituting (fTT^ . (H^Ol), (I2.22bl) . into (ITTHll and using (ITTll) shows that ah terms 
cancel, verifying the equation. 
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2.2 The Bessel kernel system 

In this case, the system of dynamical equations is given jj^ by 

dxj _ 1 {xjVk - yjXk)xk 



dak 2 ttj — Uk 

dyj ^ 1 {xjyk - yjXk)yk 
dak 2 aj — ak 



j + k, (2.23a) 

j + fc, (2.23b) 



This is again a compatible system of nonautonomous Hamiltonian equations generated 
by the Poisson commuting Hamiltonians fljC^ ■ defined in ()2.7|) . provided the Poisson 
brackets are defined by 

{xj, yk\ = —Sjk, {xo, yo} = -4. (2.24) 

aj 

There again exist two additional conserved quantities for this case. The first is defined 
by 

G^ :=xo + ^y^o+yo + ^Qi, (2.25) 

as may be seen directly by differentiating with respect to the a^-'s, using (|2.23p . Since all 
the quantities appearing in (|2.25|) vanish in the limit {aj — > 0, V j}, this difference must 
vanish, and therefore the invariant relation 

1 1 

xq = --^yl-yo- -^Qi (2.26) 

is satisfied for this solution. The second conserved quantity is 

2n 

G^ := yo + 4 ^ a^Glj = yo + ^R^i 

= yo-\{o? + xo) Qo + ^Qi + yoSo + 4Po + QqPq - Sl (2.27) 

Again, due to the Poisson commutativity of the Hamiltonians defined in ()2.7p . the Hamil- 
tonian vector fields generating the aj deformations when applied to the term R^i give 
zero, and therefore only the explicit dependence of this term upon the parameters need be 
taken into account when verifying that differentiation of the sum gives zero. Since all the 
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quantities appearing in (|2.27)) vanish in the hmit {oj — > 0, V j}, the invariant G^ must 
also vanish on this particular solution, and we therefore have the relation 



2/0 = -4i?f 1 = -4PiF, 



= \{a'' + xo) Qo - \Qi - yoSo - 4Po - QoPo + ^o- (2-28) 

The quantities R^i, Ra2 ^^^ given by 

2n 



Rl^=V^F!^ = Y,ci,G, 



B 



■4 (a' + 2;o) Qo + ^Qi + lyo^o + P^ + IQo^o - \sl (2.29a) 

lb lb 4 4 4 



i?a,2 — ^2-^a — 2^ "jG'a 



2n 



= - -^ (a^ + xo) Qi + -^Q2 + jyo^i + i^i • (2.29b) 

lb lb 4 

It again follows from the Poisson commutativity of the Hamiltonians {G^ ■} that the 
derivatives of R^i and R^2 with respect to the parameters are given by their explicit 
dependence on these parameters, and hence 

VfF^ = V,Rl, = ^Qu (2.30a) 

lb 

V2V,F^ = V2R^i = ^Q2- (2.30b) 

lb 

From (|2.23|) . application of Pi to {Qi, S'i,xo,yo} gives 

-DiQi = Qi + 4Si + lyoQi, V^Si = Si + \{a' + xo) Qi - ^Q2 + 2Pi, (2.31a) 
Pixo = -5i, Viyo = --Qi. (2.31b) 

Further application of Vi , using H2.20a|l , (|2.31|1 , and (j2.2b|l therefore gives 



VtF^ = -i^l + ^)Q, + -Su (2.32a) 

^1^" "ibV +T + T + yj^^ 



+ -Pi + { - + —]Si Qi Qa- (2.32b) 

2 V2 8 y b4^^ 32^ ^ ^ 

Substitution of (I2.29b|) . dOHl) . (ITH^ into (fTTTl) . and use of (IT^ to replace the term 
—AViF^ by 2/0) and H2.2b(l to eliminate xo, shows that all the terms cancel, verifying the 
equation. 
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3 Deduction of the r-function equations 
from isomonodromic deformations 

In this section, we show how the full hierarchies of equations derived in |21 El may be 
deduced from the Hamiltonian systems (|2.12jl . (|2.23|) and also how the corresponding 
hierarchy is deduced for the case of the sine kernel. The key step is to recast these 
systems as isomonodromic deformation equations for an associated differential operator in 
an auxiliary spectral variable z E P^, having rational coefficients with poles at the points 
{z = ttj}, and to interpret the Fredholm determinants r"^, r and t^ as isomonodromic 
T-functions. 

3.1 The Airy kernel isomonodromic system 



The Hamiltonian system (|2.12j) implies that the compatibility conditions 

§^ = lAiiM. , ^ ,, (3,., 

dak aj - ak 

^ = [a,B + C,A,]-Y^^^^^, (3.1b) 

^ = [B.^,1 (3.1c) 

are satisfied for the following overdetermined system [Sj 

^ = X^W^^ (3.2a) 

--^^*^ j = l,...2n, (3.2b) 



daj 



2n ^ 



X^(z) :=zB + C + V ^, (3.2c) 

where \I' (z, ai, . . . a^n) is a 2 x 2 matrix, invertible where defined, and 

This implies the invariance of the monodromy of the operator ^ — X (z) under changes 
in the parameters {aj}. In view of eq. 1)2. 6() . according to the constructions of '12' '13', the 
Fredholm determinant r is just the isomonodromic r-function of the system (|3.H) - (|3.2j) . 
Now define the sequence of 2 x 2 matrices 
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where the quantities Pm, Qm, Sm were defined in (|2.1H) . Expanding X (z) for large z 
gives 

oo „ 



X^iz) = zB + C+Y.^- (3-5) 

Since 

Gf=^-Tes,=a,tr[{X^)'{z)), (3.6) 



and 



G^ = ^res,=oo^tr((X^)^(z)), (3.7) 



we have 



itr((X^)^(z))=z + G^ + j:-^, (3.8) 



2 

m=0 

where 

2n , m— 1 



Ri := |] af G/ = tr {BBm+i + CB^) + ^ tr ^ B^Bm^k-i (3.9) 

j=l k=0 

(with the last term absent if m = 0) are the quantities defined in (|2.9|) . 

Using the fact that the Hamiltonian vector fields generating the aj deformations give 
zero when applied to the G^s, and hence also the R^^^ it follows that the effect of 
applying the operators V^ to R^ gives just the explicit derivatives, 

m—l 

VkRi = {m + l)tT{BB.^+k)+mtT{CB.^+k-i) + Y, I ti {Bi+k-iBrr,-i-i) (3.10) 

1=1 

(with the sum in the last term absent if ?Ti = and the second term absent if m + A; = 0). 
Applying the operator Vm to ^"^, using (|3.2b|) and (|3.4|) gives the sequence of equations 



oo „ 
fc=0 

The compatibility of these equations with 1)3. 2a() implies the following equations for the 
matrices {Bm,C}. 

m—l 

T^kBm = mBm+k-l + [C,Bjn+k] + [-B,i?m+fc+l] + / ^[Bj, Bm+k-l-l], (3.12a) 

1=0 

VkC=[B,Bk], k,m£n (3.12b) 

(where the first term of (|3.12aj) is absent if m + k = and the last term is absent if tti = 0). 
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The strategy for deriving the hierarchy of equations for r is to now choose a k-va- 
lue (fci) in p.lOf) . ()3.12|) and use these equations, together with (|3.9|) to express all the 
relevant matrix elements of the -Bm's for m < k in terms of the R^s for k < ki and 
the corresponding T>kS applied repeatedly to them. Equations ()3.12|) . for k = ki may 
then be expressed entirely in terms of these quantities, and hence in terms of repeated 
applications of the operators P^ to F = Inr . An essential step in this procedure is 
to also eliminate the additional variables xq, yo from the equations through use of the 
invariant conditions H2.14() . H2.19() . 

For example, choosing ki = 1, we note that for m = 0, eq. (|3.9|) reduces to (|2.15|) while 
for A; = 0, 1 and m = 0, (pTTn]) reduces to ^TKlji and for A; = 0, m = 0, eqs. (pTT^ give 
H2.21a() . H2.21b|) . Combining these with the invariant relations (|2.14l) . (|2.19() allows us to 
express the relevant matrix elements of C, Bq and Bi as 

XQ = VoR + i^^ yo = R, (3.13a) 

Qo = -iVoR, So = -2RVoR - V^R, (3.13b) 

Po = ^R- ^T^lR - RVlR - ]^{VoRf - R^VqR, (3.13c) 

Qi = -2R - Q{VoRf - VlR. (3.13d) 

Substituting these in eq. ()3.12b|) for k = 1 gives (|2.16|) . Similarly, eq. (|3.12a|) for /c = 1, 
m = and eq. ()3.9|) for m = 1 produce the following expressions for the relevant matrix 
elements of Bi and B2. 

Si = -ViVqR -R^ - 3R{VoRf - RV^R, (3.14a) 

Q2 = -2Ri - ViVlR - 2{VoR){ViR) - RVqR - ^RV^VoR - ^{VqR) (V^) R 

+ i (p2^)' - Ir' - ^R'^lR - 7CDoRf - Ir'CDoR?. (3.14b) 

Substitution of (|3.14bj) in eq. (|3.12a|) (or ^TTU^ ) for k = 2,m = 0, thus gives 

4ViR - 2Ri - ViVlR - 2(VoR){ViR) - RVqR - -RViVqR - -{VoR) (P^) R 

+ \ {VlRf - ^i?3 _ Ir^vIR - 7{VoRf - \R\VoRf = 0. (3.15) 

as the next equation of the hierarchy. The remaining equations may similarly be expressed 
in terms of the derivations 2?^ acting upon F . 

3.2 The Bessel kernel isomonodromic system 

The Bessel kernel case is so similar to the above that only the pertinent equations will be 
given, without repeating any details of the procedure. Define for this case, the matrices 

2n 

X^{z) -.= 3 + ^^^ + y ^^, (3.16a) 



where the Aj^s are again defined as in ()3.3a|l . 
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The Hamiltonian system (|2.2^-{)) implies that the compatibihty conditions 

J + K (3.17a) 



dAj _ [Aj,Ak] 



dak aj - Uk 



■J 

2n 



a,^—- = [Ca + ajB, Aj] - ^ _ , (3.17b) 



oaj ■f— ' aj — ak 



^ = [B,A,l (3,17c) 

are satisfied for the system 

^-=X^(z)M/^, (3.18a) 

-K— = '—^'', j = l,...2n, 3.18b 

OUj Z — ttj 

where '^^{z,ai, . . . ,a2n) is again a 2 x 2 matrix, invertible where defined. This again 
imphes the invariance of the monodromy of the operator ^ — X (z) under changes in 
the parameters {aj}. In view of eq. (|2.7)) . the Fredholm determinant t^ is again an 
isomonodromic r- function for the system (|3.17|) ~ H3.18|) . 

Defining the sequence of 2 x 2 matrices {Bm, m G N} as in (|3.4|) . and expanding X (z) 
for large z gives 

^^(-) = 5 + V + ^iSr' (3-19) 

m=l 



and 

where 

<i = ^ (C^ - Go^ - «') + ^ tr (C72 + 25i?i) , (3.21a) 

-. Tn—2 

R^^ = ir (BBm + CaBra-i) + :t tr ^ BkBm-k-i, m > 2 (3.21b) 



fc=i 



are the quantities defined in (|2.10|) and Gq , G^ are the conserved quatities defined 
in (|2.25jl . (|2.27jl . which vanish on the particular solutions defined by (|2.1j) . 

The fact that the Hamiltonian vector fields generating the aj deformations give zero 
when applied to the G^ 's, and R^„^^s again implies that the effect of applying the 
operators Vk to the -R^^'s is to evaluate only explicit derivatives with respect to the 
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parameters, giving 



VkR"^^ = mtr [BB^+k-i) + (m - 1) tr {Co,Bm+k-2) 

m—2 

+ Y,^^T^iBi+k-iBrn-i-i), m>2 (3.22) 

1=1 

(with the sum in the last term absent if m = 2). 

Applying the operator Vm to ^^, using p.4() and ()3.18b|) . again gives the sequence of 
equations 

oo „ 



P^^s = - ^ ^^Sr^"". "^ e N. (3.23) 

fc=0 

whose compatibility with (|3.18aj) implies the following equations for the matrices {Bm, Ca\, 



m—l 



T^kBm = mBm+k-1 + [Ca, B„i+k-l] + [B , -Bm+fc] + /^ [Bl, Bm+k-l-l], (3.24a) 

1=1 
VkCa = [B,Bk], k,meN, m > 1. (3.24b) 

The hierarchy of equations for t^ is derived in the same way as for the Airy case. For 
example, eqs. (|3.21() for A; = 2 reduce to 1)2. 29(1 . while (|3.22|) for A; = 1, 2, m = 1 reduces 
to (|2.3Uj) . and eqs. (|3.24|) for A; = 1, 2, m = 1 give (|2.31|) . Combining these with the 
invariant relations (|2.26l) . (|2.28|1 allows us to express the relevant matrix elements of Cq, 
Bi and B2 as 

xo = -4(Pi<i + «i)'-4iif,i), yo = -4<i, (3.25a) 

Qi = 16Pii?f^i, (3.25b) 

Si = SR^^^ViR^^ - 4Pii?f 1 + 4VfRl^, (3.25c) 
Pi = i?f,2 + a^R^^, + 4 «,)' Piiif,, - 4 (Pi<i)' 

+ ARliVlRli - P2<i, (3.25d) 

Q2 = 16P2i?f,i. (3.25e) 

Substituting these in eqs. ()3.24|) for k = 2 gives (|1.11() . Similar calculations for higher 
values of k yield the further equations of the Bessel hierachy. 

3.3 The sine kernel system 

For this case, the quantities defined in (|2.1aj) - ()2.1b|) satisfy the system of dynamical equa- 
tions defined in ^2 ^\ 

dxj _ 1 {xjyk - yjXk)xk 



duk 2 ttj — Ofc 

dyj ^ 1 {xjVk - yjXk)yk 
dak 2 Oj — Ofc 



j + k, (3.26a) 

(3.26b) 
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dxj 1 "A {xjVk - yjXk)xk ^Q Oft ^ 

a =7;}^^ + 2?/i, (3.26c) 



fc=l 



M = 1 y ^^^-^^ ~ y^'^'^y' - ^x,. (3.26d) 

This is again a compatible system of nonautonomous Hamiltonian equations, generated 
by the Poisson commuting Hamiltonians {G^} defined in ()2.5p . They imply that the 
compatibility conditions 

1^ = lAtM, , ^ ,, (3.27a) 

dak CLj - Qk 

2n 



daj " -' ■^ ttj -at 



[Bs.A,]-Y,''^2^, j^k (3.27b) 



are satisfied for the system 



dz 



X^(z)^^, (3.28a) 



-— = ^^^, j = l,...,2n, (3.28b) 

oaj z — ttj 

where 

^ A- 

X^{z):=Bs + y^^^, (3.29a) 

^^■■=[-2 f)' (^-^^b) 

with the Aj''s again defined as in 1)3. 4p . As in the previous cases, this implies the invari- 
ance of the monodromy of the operator -q^ — X {z). In view of eq. (|2.5j) . the Fredholm 
determinant t^ is an isomonodromic r-function for the system ()3.27|) - H3.28|) . 
Expanding X {z) for large z gives 

00 „ 

X'{z) = Bs + Y.^^^ (3-30) 

r?i=0 

with the matrices {Bm-, m G N} again defined as in (|3.4() . and 

-. 00 jyg 

itr((X-)^.))=-vr^+5:-^, (3.31) 

m=0 

where 

2n ^ m— 1 

^m := Yl <^J = *^ (^sBm) + 2 t^ E BkBm^k^i, m G N. (3.32) 

i=i fc=o 
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Applying the operators Vk to i?^ again just differentiates explicitly with respect to the 
parameters, giving 

m—l 

VkRi = miT{BsBm+k-i) + J^ /tr (Si+fc_iB^_i_i) (3.33) 

1=1 

(with the first term absent if /c + m = and the sum in the last term absent if ?tt, = 0) . 
Applying V^ to ^^ , using H3.30() and (|3.28b|) . gives the sequence of equations 

oo 

P^M/^ = -^-g^^^, mGN, (3.34) 

fc=0 

whose compatibility with (|3.28a() implies the following equations for the matrices {i?m}, 

m—l 
T^kBm = mB.fn+k-l + [Bs, Bm+k] + 2_^ [Bh Brn+k-l-l]- (3.35) 

The hierarchy of equations for r is derived in the same way as for the Airy and Bessel 
cases, except that we no longer have two conserved quantities like Gq ' , Goo' • To derive 
a closed system of equations, we are obliged to include two further dependent variables r^, 
which we choose as the nonvanishing entries of the matrix [Bs, Bq]t^ , 

4 ■■= (^Po - yQo) T^, r' ■■= Sot'. (3.36) 

The remaining component of Bq, which cancels in the commutator [Bs, Bq], is 

R^ = tr {BsBo) = Po + ^Qo = 0, (3.37) 

where the first equality follows from choosing ?7i = in (|3.32|) . This provides a single 
conserved quantity that vanishes for the particular solution defined by (|2.1aj) - H2.1b|) . 

To derive the hierarchy of r- function equations, we first combine eqs. H3.36() - ()3.37|) . 
which allows us to express the matrix elements of Bq as 

^0 = ~^^' ^0 = i75' ^0 = :^- (3.38) 

Eq. (|3.35|) for A; = 0, m = gives 

VoPo = -7:^So, VoQo = ASo, VqSq = 2Po - —Qo, (3.39) 

and substituting (I071) . (lOHl) in (lOIII) gives 

Vot' = 0, (3.40a) 

VoT^ = t1, VotI = -47rVf . (3.40b) 
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These equations are the lowest ones in the sine kernel hierarchy; note that they are linear 
because of the vanishing of the invariant Rq. To obtain higher, nonlinear equations, we 
first note that eq. (|I1,S2|) for m = 1 gives 

Rf = Pi + \Qi + \ {Si - QoPo) , (3.41) 

while (|3.33|) for A; = 0, 1, ?n, = 1 reduces to 

VoRf = Rq = 0, (3.42a) 

ViRf = Pi + —Qi. (3.42b) 

The first of these just gives the equation 

VqVit^ = 0, (3.43) 

which already follows from (|3.4Ua|) . The second, combined with eq. 1)3.41(1 and eq. ()3.35|1 
for A; = 1, TTi = gives the further equation 



r^Vlr^ - {V,rsy = r^V.r^ - t- (rf )^ - -^ (r|)^ . (3.44) 



,2 <;',„ c; 1 / g\2 1 / c:;\2 

4 ^'^~' ~ 167r2 
Equation ()3.35|) for fc = 1, m = gives 

PiSo = 2Pi - yQi, PiPo = -^'5i, ViQo = ASi. (3.45) 

Solving these, together with p.42b|) . gives the following expressions for the matrix entries 
of Si: 

2 ViT^ 1 /rf \^ 1 /r|\^ 1 /rf \ 

^^ = 2^ - 8 (t^J - 3^ \^) + 4^^ y^) ' (^-^^^^ 

1 _ (r^ 



Combining eq. (|3.35j) for (A; = 1, m = 1) and for (A; = 2, m = 0) gives 

P2go = ^iQi-Qi, V2Po = ViPi-Pi, V2So = ViSi-Si, (3.47) 

Substitution of (|3.38p . (|3.46p into (|3.47p gives the next equations of the hierarchy. Re- 
peating this procedure for higher (k,m) values similarly generates the higher equations. 

4 Classical i^-matrix approach and relation 
to isospectral flows 

In m 0] , a key step in deriving the hierarchies of equations for the Fredholm determi- 
nants r and t^ was to begin with certain bilinear equations satisfied by KP r-functions 
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with respect to the flow parameters {^1,^21 • • •} and to then use Vh'asoro constraints to 
replace the tm-derivations at vanishing t- values by the operators 'Dm- In this section, we 
show how the classical i^-matrix approach to the underlying isomonodromic deformation 
equations developed in |H] provides a direct link with commuting isospectral flows in the 
loop algebra s[(2), without the requirement that these arise as reduced KP flows. This fits 
into the broader framework of commutative isospectral flows in loop algebras with respect 
to the rational i?-matrix Poisson (or Adler-Kostant-Symes) structure |15[ HI IH] (and 
allows us to include the sine kernel case, which does not appear as a reduced KP flow). 

First we recall 012! that the isomonodromic deformation equations (|3.1|) . (|3.17j) . (|3.27|1 
may be viewed as Hamiltonian equations on the space of sets {^j}j=i,...,2n of s[(2) ele- 
ments, with respect to the Lie Poisson bracket, extended in the Airy and Bessel cases by 
the canonical variables (xo,yo)- (The particular form (|3.3aj) for the Aj's just represents 
a canonical parametrization on the symplectic leaves for which the Casimir invariants 
{tr74^} all vanish.) The formulae (|3.2c|) . (|3.1()a|) . (|3.29ajl deflne a Poisson embedding of 

this space into the space s((2) Jj of rational, traceless 2x2 matrices depending rationally on 
the auxiliary loop variable z, with respect to the Lie Poisson bracket on s[(2) corresponding 
to the Lie bracket: 

[X,Y]n:=l[RX,Y] + ^[X,RY], (4.1) 

where 

R:=P+- P_ (4.2) 

is the classical R-matrix, given by the difference of the projection operators 

P+ : ^1(2) ^ ;'l+(2), P+ : ^'1(2) ^ ^'(+(2), 

P^:X^X+, P_:X^X_ (4.3) 

to the subalgebras s[+(2), sl-{2) consisting respectively of the nonnegative and negative 
terms in the Laurent expanson of X(z) for large z. The space s((2)|j is identifled as 
a subspace of 5 ((2) through the trace-residue pairing 

{X, Y) := res,=oo tr {X{z)Y{z)). (4.4) 

In this setting, the isomonodromic deformation equations 1)3. 1|) . ()3.17|) . 1)3. 27(1 may all 
be expressed in the form 

- = -[(dG,)-,X] + A_^, (4.5a) 

{dG,). = -^^, (4.5b) 

z aj 

where X denotes X'^, X or X^, and Gj denotes G^, Gf or G^ j respectively. Viewing 

the Hamiltonians {Gj} as spectral invariants deflned on the space sl(2), eq. (|4.5aj) follows 
from the Adler-Kostant-Symes theorem, in view of the relations 

oaj oz 
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where ^^ denotes the derivative with respect to the explicit dependence on the parame- 
ters {aj} only. 

Rather than using the spectral invariants {Gj} as Hamiltonians, we consider the Hamil- 
tonian equations generated by the linear combinations R^, R^ or i?^„ defined in (|2.8() . 
((T^ . ((TTU)) . which are ah of the form 

V^X = -[{dR^)^,X] + ^(^^, (4.7) 

with the respective identifications for X and {Rm}- These are just equations (|3.12() . (|3.24|) 
or (|3.35|) . depending on the identification, since 

Rm = lres,=ooz"'trX\z), (4.8) 



(4.9) 



and therefore dRm, viewed as an element of s 1(2 


l), is just 


oo „ 
k=0 




implying 




fc=0 





(4.10) 

If, instead of the nonautonomous systems occurring here because of the identifications 
of the Oj's as multi-time parameters, we consider the autonomous systems generated by 
the same set of Hamiltonians {Rq,Ri, . . .}, denoting the corresponding flow parameters 
{to,ti, . . .}, the resulting equations have the isospectral form 

dX 

-— = ±[(di?„)±,X], (4.11) 

UtfYi 

where either of the projections {dRm)± may be used, since the differential dRm^ given 
by (|4.1UI) . commutes with X. Although these systems are generated by the same Hamilto- 
nians as the nonautonomous systems (|4.7|1 . they of course do not generate isomonodromic 
deformations of the operator -4- — ^{z), and in fact are not even compatible with the 
systems ()4.7() : however, they are compatible amongst themselves, generating commuting 
isospectral Hamiltonian flows. The close relationship between the autonomous and as- 
sociated nonautonomous systems implies a correspondence between the structure of the 
resulting hierarchies. 

To see this, we substitute the expressions H3.2c|) . (|3.16a|) and ()3.29a|) for X{z) and ()4.1U|) 
for dRrn into (|4.11|1 to obtain the systems 

dB ^~^ 

— — = [C, Bm+k] + [B, Bm+k+i] + 2_^ [Bi,Bm+k-i-i], (4.12a) 

"^ ;=o 

-- = [B,Bk], k,meN (4.12b) 

Otm 
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for X = X^, 



dBm 

dim 



m—l 



[C, -Bm+fc-l] + [B, i^m+fc] + 2^ [Bl,Bm+k-l-l], (4.13a) 



1=1 



dt 
for X = X^ and 

dB„ 



[B,Bk], k,men, m>l. (4.13b) 



m— 1 



dtr, 



[Bs, B^+k] + Yl [^'' B^m+k-i-i] (4.14) 



for X = X^ . These only differ from the equations 1)3.121) . 1)3. 24() and 1)3. 35() by the absence 
of the term mBm+k-i in the right hand side of H4.12b|) . H4.13b|) . 1)4. 14|) and the replacement 

Pm - ^ (4.15) 

for the derivation on the left hand side. The procedure for deriving hierarchies for such 
systems is well known in the isospectral context (see, e.g. [7] for details); the recursive 
procedure used in Section 3 above is just the analog of this approach applied to the 
isomonodromic systems ()3.12|) . ()3.24|) and 1)3. 35(1 . 

As a final point, it should be noted that almost nothing in the derivation of the r- 
function equations of Sections 2 and 3 depended on the fact that the specific r-functions 
involved were equal to the Fredholm determinants ()1.2|) . H1.4() . (|1.6() . Everything just 
followed from the general form of the isomonodromic deformation equations ()3.1j) . ()3.17)1 
and 1)3. 27() . the only features specific to the identifications of r , r^, t^ as Fredholm 
determinants being the fact that the matrix residues Aj were of rank 1 (as seen from 
the parametrization 1)3. 3a|) ) and the invariants Gq, G^, Gq , G^ vanished. By allowing 
these invariants, as well as the constants {detvlj}, to take arbitrary values, an identical 
procedure leads to equations for the r-functions of the general isomonodromic systems, 
which only differ from the ones derived in Sections 2 and 3 by the nonzero constant values 
of the two additional invariants Gq , G^ or Gq and G^. For example, eq. (|2.1H|) is replaced 
in the general case by 

vIr - AViR + 2R + A{gl^- go) VqR - 2g^ [vIr + 2RVqR) + 6{VoRf = 0, (4.16) 

where go, goo are the values taken by the invariants Gq , G^, respectively. The other 
equations of these hierarchies may similarly expressed in a way that allows arbitrary 
values for these constants. 

Acknowledgements 

The author would like to thank P. van Moerbeke, C. Tracy and H. Widom for helpful 
discussions relating to this work. Research supported in part by the Natural Sciences and 
Engineering Research Council of Canada and by the Fonds FCAR du Quebec. 



Equations for Fredholm Determinants Appearing in Random Matrices 549 

References 

[1] Adams M R, Harnad J and Previato E, Isospectral Hamiltonian Flows in Finite and Infinite 
Dimensions I. Generalised Moser Systems and Moment Maps into Loop Algebras, Commun. 
Math. Phys. 117 (1988), 451-500. 

[2] Adler M, Shiota T and van Moerbeke P, Random Matrices, Vertex Operators and the Virasoro 
Algebra, Phys. Lett. A208 (1995), 67-78. 

[3] Adler M, Shiota T and van Moerbeke P, Random Matrices, Virasoro Algebras and Noncom- 
mutative KP, Duke Math. J. 94 (1998), 379-431. 

[4] Adler M and van Moerbeke P, Completely Integrable Systems, Euclidean Lie Algebras, and 
Curves, Adv. Math. 38 (1980), 267-317; Linearization of Hamiltonian Systems, Jacobi Varie- 
ties and Representation Theory, Adv. Math. 38 (1980), 318-379. 

[5] Borodin A and Deift P, Fredholm Determinants, Jimbo-Miwa-Ueno Tau-Functions and Rep- 
resentation Theory, math- ph/01 1 1007 

[6] Forrester P J, The Spectrum Edge of Random Matrix Ensembles, Nucl. Phys. B402 (1993), 
709-728. 

[7] Flaschka H, Newell A C and Ratiu T, Kac-Moody Lie Algebras and Soliton Equations, 
Physica D9 (1980), 300-323. 

[8] Harnad J, Dual Isomonodromic Deformations and Moment Maps into Loop Algebras, Com- 
mun. Math. Phys. 166 (1994), 337-365. 

[9] Harnad J, Tracy C A and Widom H, Hamiltonian Structure of Equations Appearing in Ran- 
dom Matrices, in Low Dimensional Topology and Quantum Field Theory, Editor: Osborn H, 
Plenum, New York, 1993, 231-245. 

[10] Harnad J and Its A R, Integrable Fredholm Operators and Dual Isomonodromic Deformations, 
Commun. Math. Phys. 226 (2002), 497-530. 

[11] Jimbo M, Miwa T, Mori Y and Sato M, Density Matrix of an Impenetrable Bose Gas and 
the Fifth Painleve Transcendent, Physica Dl (1980), 80-158. 

[12] Jimbo M, Miwa T and Ueno K, Monodromy Preserving Deformation of Linear Ordinary 
Differential Equations with Rational Coeefficients I, Physica D2 (1981), 306-352. 

[13] Jimbo M and Miwa T, Monodromy Preserving Deformation of Linear Ordinary Differential 
Equations with Rational Coeefficients II, HI, Physica D2 (1981), 407-448; D4 (1981), 26-46. 

[14] Mehta M L, Random Matrices, 2nd edition. Academic, San Diego, 1991. 

[15] Olshanetsky M M, Perelomov A M, Reyman A G and Semenov-Tian-Shansky M A, Integrable 
Systems II, in Encyclopaedia of Mathematical Sciences, Vol. 16, Springer- Verlag, Berlin, 
Heidelberg, New York, 1994. 

[16] Palmer J, Deformation Analysis of Matrix Models, Physica D8 (1994), 166-185. 

[17] Tracy C A and Widom H, Introduction to Random Matrices, in Geometric and Quan- 
tum Methods in Integrable Systems, Springer Lecture Notes in Physics, Vol. 424, Editor: 
Helminck G F, Springer-Verlag, N.Y., Heidelberg, 1993, 103-130. 



550 J Harnad 

[18] Tracy C A and Widom H, Level Spacing Distributions and the Airy Kernel, Commun. Math. 
Phys. 159 (1994), 151-174. 

[19] Tracy C A and Widom H, Level Spacing Distributions and the Bessel Kernel, Commun. Math. 
Phys. 161 (1994), 289-309. 



